EXISTENCE OF SOLITARY WAVES IN DIPOLAR QUANTUM 

GASES 



PAOLO ANTONELLI AND CHRISTOF SPARBER 

Abstract. We study a nonlinear Schrodinger equation arising in the mean 
field description of dipolar quantum gases. Under the assumption of suffi- 
ciently strong dipolar interactions, the existence of standing waves, and hence 
solitons, is proved together with some of their properties. This gives a rig- 
orous argument for the possible existence of solitary waves in Bose-Einstein 
condensates, which originate solely due to the dipolar interaction between the 
particles. 



1. Introduction 

The experimental realization of Bose-Einstein condensation (BEC) in dilute 
gases in 1995 [3] has marked the beginning of a new era in atomic physics and 
quantum optics. Ever since then, continuous efforts have been undertaken to ex- 
tend BEC physics towards new regimes offering different theoretical and experi- 
mental challenges. A particularly interesting research field concerns the study of 
solitary waves within BECs, see e.g. [1] and the references given therein. More 
recently, so-called dipolar BECs, i.e. condensates made out of particles possessing 
a permanent electric or magnetic dipole moment [24] , have received much attention. 
This is due to the fact that the additional dipolar interactions between particles are 
both long-range and non-isotropic and therefore crucially influence the ground state 
properties, stability, and dynamics of the condensate, see [16] for a broad review 
of this subject. In addition, the possibility of a novel class of solitary waves within 
such systems has been discussed in e.g. [13, 21, 23]. Motivated by these reports, it 
is the aim of this paper to rigorously prove the existence of solitary waves within 
dipolar quantum gases. 

To this end, we shall be concerned with the mean-field description of (dilute) 
dipolar quantum gases, based on an (augmented) Gross- Pitaevskii equation, cf. 
[16] for the validity of such a description. Following [27, 28], we shall describe the 
time-evolution of such systems by the following Gross-Pitaevskii type model 

(1.1) ihdti} = /^ip + g\il)\'^ilj + a'^{K i, e]R,xeR^ 

2m 

where g = 4:Trh?Na/m, for a G N number of particles, with mass m > and 
scattering length a G M. Finally cr^ ^ denotes the strength of the dipole moment. 
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In (1.1) we also denote by " * " the convolution w.r.t. x between the loeal density 
p = IV'P and 

N 1 — 3cos^ 9 
(1.2) K{x) = 

where 9 = 9{x) is the angle between x € and a given (fixed) dipole axis n e M^, 
such that |n| = 1, i.e. 

x-n 

cos 6 = , , . 



\x\ 

The non-local potential K*p describes long-range dipolar interactions between the 
particles, whereas the local (cubic) nonlincarity oc \'ip\'^'ip describes the usual contact 
interaction between particles (which is short-range, isotropic and characterized by 
the scattering length a e M). 

For the upcoming mathematical analysis it will be more convenient to rescale 
(1.1) into the following dimensionless form 

(1.3) idt^ = —Ai, + AilVI'V- + HK * |Vl')V', 

where Ai, A2 £ M are some constants (depending on the physical parameters given 
above), which describe the strength of the two nonlinearities, respectively. The 
time-dependent equation (1.1), in the presence of an additional (quadratic) confin- 
ing potential, has been rigorously analyzed by Carles, Markowich and the second 
author in [8]. Several existence and uniqueness results are discussed in [8], as is 
the possible occurrence of finite time blow-up of solutions, which physically corre- 
sponding to the possible collapse of the BEC [20], see Remark 1.3 below for more 
details. 

Prom the mathematical point of view, it is well known (see e.g. [9]) that the 
possibility of finite time blow-up is strongly linked to the existence of standing wave 
solutions to (1.3), i.e. 

ip{t,x) = e'"*u{x), LoeR. 

Obviously, for tp{t,x) to be a solution of (1.3), the profile u € L^(R^) has to solve 
the following nonlinear elliptic equation 

(1.4) MH'^u + X2{K * \u\'^)u + uju = 0, 

which will be the main object of study in our work. To this end, we first recall that 
the total energy associated to associated to (1.4), is given by 

(1.5) E{u):= [ )-\Vu\''+'^\u\^+'^{K*\u\^)\u\''dx = T{u) + V{u), 
where 

(1.6) T{u) ~ \ j \Vufdx, 
is the kinetic energy, and 

(1.7) V{u):= j ^\u\' + ^{K*\u\^)\u\^dx, 

is the nonlinear potential energy. At this point it might be tempting to study (1.4) 
via minimization of the energy E(u). However, it is well known, that even without 
the dipole nonlincarity i.e. A2 = 0, this approach fails, since, on the one hand, the 
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energy functional in three spatial dimensions is found to be unbounded from below 
in the case Ai < (by invoking the Gagliardo-Nirenberg inequality). On the other 
hand, when A2 = 0, Ai > the minimization problem becomes trivial. In other 
words, equation (1.4) turns out to be L"^ -supercritical in the sense of [9] and this 
problem is in fact enhanced by the presence of the dipole nonlinearity. Nevertheless, 
we shall still follow a variational approach for studying the existence of solutions 
to (1.3). This approach is based on the chocie a suitable functional (see Section 
3 below) which has first been introduced in [26] in the case of local nonlinearities. 
Since then, it has proved its use in different applications, in particular for water wave 
models of Davey-Stewartson type [10, 11, 22], which have a similar mathematical 
structure as the model we consider here. 



Theorem 1.1. Let Ai, A2 € M be such that the following condition holds: 
(1.8) Ai < < 



47r 

yA2, «/A2>0, 



-g-A2, i/A2<0. 

Then there exists u G if^(R^) solution to (1.4) with corresponding w > 0. Further- 
more the solution u satisfies the following properties: 

(1) u is real-valued and u{x) ^ 0, Va; G M^. 

(2) u is radially symmetric in the X\,X2— plane and axially symmetric with 
respect to the x^ — axis, i.e. 

u{xi cosa + X2 sin a, —xi sina + X2 cos a, X3) = u{xi,X2, X3), 

V a G [0,27r], andu{xi,X2, —X3) = u{xi,X2,xs). 

(3) The energy of u satisfies E{u) = ^T{u) > 0. 

(4) u e H'(R^) for all s^l. 

(5) There exist positive constants Ci,C2, such that: 

e^il"^l(|M(.x)| + |V?i(.T)|) s; C2, Va; G 

The assumptions on Ai , A2 G IR. can be interpreted as the necessity of a sufficiently 
strong dipolar nonlinearity. Note that the existence of steady states is guaranteed 
even in situations where Ai > 0, i.e. in the case of a repulsive (defocussing) cubic 
nonlinearity. The appearance of steady states in this regime is therefore solely 
due to the presence of the dipolar interaction and can not be reproduced in a 
conventional BEG. Moreover, by invoking the Galileian-symmetries of (1.3), the 
existence of steady states directly implies the existence of dipolar solitons in the 
form 

(1.9) ip{t, x) = u{x + K^)e*'^*e-'''(^+'^*)/^ k G R. 

where u is the smooth, exponentially decaying (as |a;| — > 00) non-negative profile 
guaranteed by Theorem 1.1. In the usual language of solitary waves, such a solution 
to (1.3) is considered to be a bright soliton [1]. Our work therefore provides a 
rigorous mathematical basis for the existence of dipolar solitons, as studied in [13, 
21, 23] (see also [16] for a broader discussion and [6] for a closely related physical 
system). 

The main drawback of Theorem 1.1 is that it leaves uniqueness as an open 
question. In the usual case of a single cubic nonlinearity, uniqueness of (positive) 
solutions is strongly interwoven with the fact that u = u{\x\) is found to be radially 
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symmetric [12] (the proof is based on symmetric re- arrangement techniques). This 
radial symmetry is broken in om- case due to the influence of the dipolar nonlinearity 
and hence we can not conclude uniqueness. 

Remark 1.2. In physical experiments one needs to confine the gas by electromag- 
netic traps in order to achieve sufficiently low temperature. In order to mathemat- 
ically describe the trapping potential, the Gross-Pitaevskii equation usually carries 
an additional term, namely 



where the V^xt is assumed to be of the following form Vc^t (x) — , i.e. a harmonic 
oscillator confinement. Obviously, the presence of 14xt accounts for the existence 
of steady states even in the linear case Ai = A2 = (where the Hamiltonian 
admits countable many eigenvalues). The situations with trapping potential there- 
fore has to be clearly distinguished from the one considered in Theorem 1.1 above. 
In particular, if one assumes the presence of a confinement and in addition, say, 
Ai ^ ^ 0, one can easily obtain the existence of stationary states in trapped 
dipolar BEG by following the arguments given in e.g. [19] (see also [5]). That is, 
by minimizing the corresponding energy functional E(u), which is now obviously 
bounded from below since all terms within E{u) are positive (see also Remark 2.3 
below) . Clearly, the presence of T^xt breaks the Galileian-synimctry of the model 
and thus solitary waves of the same kind as given by (1.9), can only be observed 
in an actual physical experiment, when the trapping potential is turned off and the 
BEG is allowed to evolve only under the influence of nonlinear effects. 

Remark 1.3. Under Assumption (1.8) on Ai,A2, finite time blow-up (and hence 
collapse of the condensate [20]) can occur for the time-dependent equation (1.3). 
Indeed by a straightforward calculation (see also [8, 20]) one obtains the virial 
identity 



Invoking the classical argument of Glassey [14] yields blow-up of solutions to (1.3) 
in finite time, provided (1.8) holds true and the initial energy is negative (see [8] 
for a possible construction of such initial data). Note, however, that the solitary 
wave solutions constructed above correspond to initial data with positive energy 
(see assertion (3) of Theorem 1.1). They do not blow-up in finite time but most 
probably are unstable w.r.t. small perturbations (an issue that can be overcome in 
experiments by creating effective lower dimensional systems, cf. [4]). 

The paper is now organized as follows: In Section 2 below we shall as a flrst step 

derive necessary conditions for the existence of standing waves. That these waves 
in fact exist is then proved in Section 3 and we conclude the work by giving the 
remaining details for the proof of our main theorem in Section 4. 
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2. Necessary conditions for, existence 

In this section we shall prove several conditions which are necessary for the 
existence of solutions to (1.3). Note that without loss of generality, we can impose 
n = (0, 0, 1). In this case K{x) becomes 

(2.1) Kix) = 

and we shall assume K to be of this form from now on. In the following we also 
denote the Fourier transform of a function f{x) by 



JR3 



We can then recall Lemma 2.1 of [8] concerning basic properties of the nonlocal 
potential K * \u\^. 

Lemma 2.1. The operator IC : f t-^ K * f can be extended as a continuous operator 
on LP(M.^), for each 1 < p < oo. Moreover, the Fourier transform of K is given by 



(2.2) ^(0 = ^^(3cos^e-l) = ^^'^^'-^^'-^' 



3 ^ ' 3 V I^P 

where 6 = 0(^) denotes the angle between ^ G and the dipole axis n = (0, 0, 1). 

Formula (2.2) implies that K G L°°(M^) and thus K * f clearly defines a con- 
tinuous operator K, : L'^(JS?) L^(R^). In the following lemma we shall derive 
two Pohozaev-type identities, which have to be a-priori satisfied by any solution to 
(1.4). 

Lemma 2.2. Let u e iJ^(M^) be a solution to (1.4). Then the following identities 
hold 



(2.3) / \\Wu\^dx = 2,u [ \u\'^dx 



(2.4) / ^\uf + ^{K*\u\^)\u\^dx = -2uj I \u\''dx. 

Jm.s 2 2 



Moreover we also have that V{u) = — |T(u) = —2E{u). 

Proof We first multiply equation (1.4) by x ■ Vu and then integrate by parts. 

Straightforward calculations yield 

(2.5) = j^^ -\\SJu\^ - l^ut - \\2{K * \u\^)\u\^ - luj\u\^dx. 
On the other hand, multiplying equation (1.4) by u, we obtain 

(2.6) = / \\Vu\^ + M\u\^ + \2{K*\u\^)\u\'' + Lo\u\^dx. 

Jm? 2 

By combining the identities (2.5) and (2.6), we infer 

0= / ~\ut -^{K *\u\^)\u\^ -2u:\u\Hx, 

which is nothing but (2.4). In view of (2.6), this also yields (2.3). Finally, the 
definitions of the kinetic and potential energy given in Section 1, together with 
(2.3) and (2.4) directly imply V{u) = -|T(u) = -2E{u). □ 
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From the Pohozacv-type identities (2.3), (2.4) we can derive the following nec- 
essary conditions for the existence of solutions to (1.4). First of all, (2.3) obviously 
requires co > 0. Consequently, (2.4) implies 

(2.7) / ^\u\* + ^(K*\uf)\u\'dx<0. 

Denoting by p(^) = |wP(^) the Fourier transform of jup, Plancherel's theorem 
allows us to rewrite the left hand side of this inequality in the following form 

/ Xi\u\^ + X2{K*\u\^)\u\^dx= [ Xip\0 + K{0p\0d^- 

In view of formula (2.2) we note that in fact K{^) G [—^, In the case of 
A2 > 0, this implies 

which consequently requires Ai < 42^A2, as a necessary condition for the existence 
of solutions, cf. the first line of (1.8). Similarly, in the case when A2 < 0, we obtain 



Ai + ^A2 ) nm, 



which yields the second condition given in (1.8), i.e. Ai < — ^A2. 

Remark 2.3. The conditions on Ai, A2 can also be understood as follows: Straight- 
forward calculations show (cf. [16]) that, for K{x) given by (2.1), the dipolar 
nonlinearity can be rewritten as 

where $ solves the Poisson equation — A$ = \u\^. The total nonlinear potential 
energy (1.7) is therefore given by 

(2.8) V{u) = \( (Ai - ^X^MUx + ^ / \d^,V^^dx, 

i.e. a term which stems from a (combined) cubic nonlinearity and a sub-critical 
term, which stems from Poisson's equation. The assumptions on Ai , A2 consequently 
ensure that the total potential energy is essentially the same as in the case of 
an attractive cubic nonlinearity. Note that (2.8) implies that in the defocusing 
situation, i.e. Ai ^ ^A2 ^ 0, the nonlinear potential energy is indeed convex (see 
also [5]). 

3. A VARIATIONAL FORMULATION 

We shall now formulate a variational problem, which will be used to ensure the 
existence of solutions to (1.4). To this end, we introduce 



(3.1) 
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where (•,•) denotes the scalar product in L'^{R^). This functional is well-defined 
for each v G H^(M.^) in view of Lemma 2.1. Moreover, J{v) satisfies the following 
scaling properties: Let 

(3.2) Vg^s{x) := qv{sx), q,s> 0, 

which implies 

Then J{v) is found to be invariant under this scaling above, i.e. J{vs,q) = J{v)- 
To this end, one checks that the nonlocal term J {K * \ v\'^)\v\'^ dx scales in the same 
way as the L*(IR^) norm, since 

/ iK*\v,,s\^)\Vg,s\^dx = q^s-^ [ K{Of (-) 

= q's-^ I KisOp^iOdC^q^s-^ f iK*\vnv\'dx, 



where p denotes the Fourier transform of p := 
the fact that -ft'(s^) = K{(^), for each s > 0, which is easily seen from (2.2). 

Lemma 3.1. Let v € H^{M.^) be a critical point of the functional J, with J{v) = 
a G M. Then v is a weak solution to 

(3.3) -l3Av + 4:a{Xi\v\'^v + X2K:{\v\'^)v) +UJV = 0, 

where uj = ||Vt;|||_2||^^||^2^ and P = 3||w||i,2||Vt;||i,2. 

Proof. Denote (3i = \\v\\i^2 and /32 = ||V?;||/^2. All we need to do is to find the 
points V G H^{M.^), where the first variation of J satisfies 6J{v) = 0. To this 
end, we calculate the variation of the nonlocal term (/C(|vp), in the form 
{2vlC'{\v\'^),r]), \fr] G C^(R^), where JC'(f) denotes the Frechet derivative of the 
operator JC{f). In order to compute /C'(/) we use Plancherel's theorem, to write 

m + V),f + V)- mf), f) = f K {2Re{frj) + |r?H dx. 
Using (2.2), we obtain 

\mf+v)j-v) - mf),f) - mf),v)\ < yii^iii^ 

and we consequently infer that the Frechet derivative of K{f) is given by (/C'(/), r]) = 
(2/C(/),j7). We therefore find that SJ{v) is given by 

3\\v\\l4'^v\\lH-Av,v) + \\Vv\\l4v\\-J{v,v) 



5J{v)= - 
+ 



iXi\\v\\l,+X2{IC{\vn\v\^)) 
WVvWUHl^^ {^X,{\v\^v,r^)+4X2{IC{\vnv,ri)) 
{Xi\\v\\l. + X2{IC{\v\^),\v\')y 



{X4v\\i.+X2{ici\vnm 

- 4a{X,{v,v) + X2{)C{\v\^v),v))). 

Since this has to hold for any r] G Cq°{M.^), we conclude 

3/3i/32At; - iS^^p^v - Aa {Xi\vfv + X2}C{\vfv)) = 0, 
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which proves the assertion of the lemma. □ 

Consequently the problem of existence of solutions to (1.4) reduces to the task 
of finding critical points, say minima, of the functional J{v). 

Proposition 3.2. Let Ai,A2 satisfy assumption (1.8). Then there exists a mini- 
mizer t;* e i?^(IR^) for the functional defined in (3.1), i.e. 

(3.4) J{v^) = j := inf J{v). 

Proof. First we note that condition (1.8) ensures that J{v) is non-negative. Hence, 
there exists a minimizing sequence {iiralnGN C iJ^(M^), i.e. 

j = lim J{vn). 

n—>oo 

Furthermore, we can rescale the sequence {(VijueN, such that 

WvtiWl^ = 1, ||Vt;„||i,2 = 1 

and thus 

(3.5) j = lim J{vn) = - ( lim Vivn))'' 

n—>oo \n—>oo / 

Since {wrijnGN is uniformly bounded in H^{M.^), we know that, up to extraction of 
a sub-sequence, there exists a weak limit in H^{R^), i.e. 

u„^w, iniJ^(]R3). 

Moreover, by the lower semicontinuity of the norm, we have 

||w,||i2 =6i < 1, \\Vv4l2 =b2 

First of all, by arguing as in [22] wo can show that both bi and 62 are strictly 
positive: Indeed, from the boundcdness of J(w„) < 00, the fact that ||w„||£,2 = 
llViirallia = 1, and assumption (1.8), we infer that the L**— norm of the sequence 
(i'n)raGN IS Uniformly bounded away from zero, i.e. Hi^nllL-i > C > 0. Hence by 
arguing as in Section 5 of [22] (which in itself is based on a result by Lieb [17]), we 
conclude that 61 ^ and 62 7^ 0. 

Thus if we can prove that {w„}„gN indeed converges (not only weakly but) 
strongly towards i;* € H^{M.'^), then clearly is a minimizer for the functional 
(3.1). To this end, it suffices to prove that indeed 61 = 62 = 1- 

We consequently define the difference w„ := w„ — for which we have 

I ||2 n— >oo ^ ,2 iirv Il2 n— >oo ^ ,2 

Hence 

j < lim JK = ^ y out \ ■ 

We now want to study the limit of 2V{wn), as n — >■ oo, which in view of (3.5) can 
be written as 

lim 2Viwn)= lim (Ai(||^i;„||l4 - ||^;„||l4) + A2((/C(^i;^),«;2) _ (/C(„2)^^2^-)^ _ 1^ 

n— foo n— >-oo J 

We first consider the term ||w„||^4 — 11^^11 ^4- Using a classical result by Brezis and 
Lieb [7] we know that 



lim (||u;„||l,4 - ||wn||l,4) 



l|4 
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In a second step we rewrite the term (/C(w^j), w^) in the fohowing way 

(/C(w;^),w^) = / wllC{vl)dx+ wllC{vl)dx - 2 VnV^IC{wl)dx 

Jr3 Jts^3 Jt^3 

= {K{vl),vl) + I wi}C{vl)dx+ I (vl-2v^Vn)lC{vl)dx 



2 / VnV^IC{w'^)da 



By Sobolev embeddings, the weak convergence w„ ^ in 
in L^'(IR^), for 1 < p ^ 3. Hence we immediately get 



implies ^ 



I 



w„lC{v^)dx > 0, 



since |l/C(w^)||^2 ^ %"lk*l|L2. Next, we want to show that 



(3.6) 
Indeed, 



2v»Vn)JC{v'^)dx 



)dx. 



R3 



{vl-2vnV»)JC{vDdx+ I vlK.{vl)dx 



K3 



K3 

+ 2 



and thus, by the Cauchy-Schwarz inequality 



/ {vl- 



2vnV^)lC{vl)dx + I vlK{vl)dx 



1/2 



R3 



wlvldx 



1/2 



w„{vn + v^)K.{vl)dx 
I WnV*lC{vl)dx , 

wlK?{vl)dx\ ' X 



/C {v^dx 



1/2 



Having in mind, that /C (cf. Lemma 2.1) is bounded in _L^(M^), we can again use 
the convergence of to^ — ^ in i^(K'^), to obtain (3.6). By similar arguments one 
can show that the term —2 J t;„t;*/C(w^)rfa; converges to zero. 

In summary, we conclude that 



- lim 2V{w^) = 2V{v.) + -, 

n—^oo J 



which consequently implies 



j < J{Wn) 



~^ 2V{v.) + \ 



Rearranging this inequality and using the fact that — j2y(w») = 6162 we obtain 
(3.7) 



1 ^ (l-6f)V2(l_;,2/3^3/2 ^^j,3 

J J 3 



On the other hand, by simple algebra we obtain 

{i-bir/\i-bi/'f/'+hbUi, 
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and thus, also the reverse inequality corresponding to (3.7) holds. Thus 

j j 3 

This identity holds if and only if, either: &i = 62 = 0, which is not allowed in our 
study, or else, if 61 = 62 = 1 • We therefore conclude that 

lk*llL2 = l, ||Vt;*||i2 = 1, 

which consequently implies the strong convergence of the minimizing sequence. In 
summary, this shows the infimum of the functional is attained at w*. □ 

4. Proof of Theorem 1.1 

Recall that the functional J{v) is invariant under the scaling (3.2). Thus we can 
choose the minimizer to be such that ||t^*||L2 — 11Vv*||l2 = /32. Hence, by 
combining the results stated in Proposition 3.2 and Lemma 3.1, we conclude that 
w* is a solution to 

(4.1) -3/3i/32A^;, + 4j {Xi\v^\^v^ + X2JC{\v^\^)v^) +u;v^ = 0, 

with u) = (3^^(32- Since /3i,/32 > are arbitrary, we can choose these parameters in 
the following way: Pi = ^ (f )"^^^' '^2 = (t)^^^) so that (4.1) becomes 

(4.2) ~Av^+4:j{Xi\v^\'^v^ + X2lC{\v^f)v^) +iov^ =0, 

where j = J(w*), as in Proposition 3.2. Hence, by rescaling u{x) — (4j)^/^u,(x), 
the existence of a weak solution u G H^{R^) to the original problem (1.4) is proved. 

It remains to prove Assertion (1) - (5) of Theorem 1.1: In order to show Assertion 
(1) we note that, by straightforward calculations 

\\V\v\\\l2^\\Vv\\l2, "^vGH'im^). 

Thus J{\v\) ^ J{v) and hence J(u,) = J(|?;*|), which implies that the minimizer 
satisfies v*{x) > 0. To prove Assertion (2) we can use Steiner symmetrization (cf. 
[25], proof of Theorem 3, and [2]. See also [18] for the a closely related topic of 
re-arrangement inequalities.). Let u"^ be the Steiner symmetrization of u around a 
plane in M^. It is straightforward, cf. [25, Theorem 3] to prove that J{w^) < J(w), 
provided that ii"* = K, where K"^ is the Steiner symmetrization of K. This means 
that we can find a minimizer u of J which has the same symmetries of K , hence it 
is radially symmetric in the (xi, 0:2)— plane and axially symmetric with respect to 
the 0:3— axis. Assertion (3) then directly follows from the identity 

E{u) = ^Tiu) > 0, 

as proved in Lemma 2.2. Assertion (4) follows from the simple observation that for 
u e H^{R'^) the right hand side of 

— ^Au = — Ai|u|^u — X2{K * |wp)t< — cou, 

is bounded in L^(M^), due to the properties of /C stated in Lemma 2.1 and the 
Gagliardo-Nirenberg inequality ]]u]]^4 ^ ||m||l2 ]] Vm|||^2- Thus we conclude that in 
fact u e if^(IR^) and an induction argument shows u e H^{R^) for all s ^ 1. 
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Finally, in order to prove the asymptotic decay as \x\ — > oo, all one has to do 
is follow the arguments given in Step 6 of the proof of [10, Theorem 2.4] (which in 
itself follows from [9] ) . 

Remark 4.1. As a by-product of our analysis we obtain, that there exists a C > 0, 
such that for any / £ H^(M.^) the following inequality holds 

-Aill/lll. - X2{K * l/p, ^ C||V/||i.||/|U., 

where the optimal constant C = C* is given by 

^ _ --^ill^*lll4 - X2{K*\v^\'^, \v^\'^) 

||VW*|||2||W*||I,2 

with t;* being the minimizer of J{v), as guaranteed by Proposition 3.2. 
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